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FOREWORD 


This final report presents work which was conducted for 
Langley Research Center (LaRC) in response to requirements of 
Contract NAS1-15819. The work presented was performed at 
REMTECH's Huntsville office and is entitled "Wake Flowfields for 
Jovian Probe." 

The NASA technical coordination for this study was provided 
by Dr. James Moss of the Aerothermodynamics Branch of the Space 
Systems Division. 
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NOMENCLATURE 


a Term defined in Equation 11; also, sonic 

velocity 

A Area 

AF Air Force 

b Intercept term in far wake shape - see Equation 9 

C Elemental concentration 

c Specific heat at constant pressure 

CEC Chemical Equilibrium Composition 

CESL Chemical Equilibrium Shear Layer 

d Diameter of probe 

F Fraction of internal species mixed with external 

species (similarity variable) 

H Enthalpy 

HYVIS Hypersonic viscous shock layer program 
m Exponent term in far wake shape - see Equation 9 

M Mach number 

M Molecular weight 

MOC Method of Characteristics 

n Term defining two-dimensional or axisymmetric 

flow 

P Pressure 


q Heating rate 

r Radial distance fro® axis of symmetry 

R Gas constant; also, Radius of probe 
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T 

u, U 
X 

y 

Y 

Y 
6 

6 * 

AY 

0 

y 
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Reynolds number based on freestream conditions 
and body diameter 

Temperature 

Velocity 

Axial distance from aft end of probe shoulder 
Radial distance from axis of symmetry 
Normal coordinate (see Fig. 3) 

Ratio of specific heats 

Boundary layer thickness 

Boundary layer displacement thickness 

Width of mixing region 

Angle, momentum thickness 

Viscosity 

Density 

Stream function defined in Equation 12 


Subscripts 


b Base 

c Cone 

d Diameter of probe 

e Edge of boundary layer 

N Neck 

r Reference 

s Shock, stagnation 
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Wake, wall 

Freestream conditions 
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Section 1 
INTRODUCTION 


The next planetary mission that will encounter aerodynamic 
heating is Galileo. After more than a 3-year journey to Jupiter, 
the Galileo probe will enter Jupiter's atmosphere with a rela- 
tive entry velocity of 48 km/sec, producing the most severe ther- 
mal environment ever encountered. In order for the probe to sur- 
vive the severe aerothermal entry conditions a large portion of 
the probe's weight must be dedicated to the thermal protection 
system. Design of the thermal protection system must rely on an- 
alytical predictions since the entry thermal environment cannot 
be duplicated in any existing ground test facility. 

The thrust of the current work is focused on the probe's 
wake flowfield. The wake produced by the probe is highly ener- 
getic yielding both convective and radiative heating inputs to 
the base of the probe. In order to calculate the radiative heat- 
ing to the base region of the Jovian probe, the entire wake flow- 
field must be defined. This study addresses the development, im- 
provement, and use of engineering tools for calculating the 
probe's wake flowfield. 

The current study is an extension of previous work (Ref. 1) 
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and includes several calculation improvements as well as develop- 
ment of flowfields for three cases. The improvements consisted 
of developing a shoulder expansion code, calculating the near 
wake recirculation zone pressure distribution and improving the 
recirculation zone model. Flowfields for three entry conditions 
were calculated and include temperature, pressure and elemental 
species concentrations. These flowfields extend radially for 2*5 
body radii and axially for 7.0 body radii. Section 2 of this re- 
port describes the math model used in the calculations. Section 
3 provides a description of the results, and this is followed by 
conclusions and recommendations in Section 4. 
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Section 2 

WAKE FLOWFIELD MODEL 

The objectives of this study required the improvement and 
implementation of an engineering model for the Jovian probe wake 
flowfield. The basis of the model was developed by Engel in Ref. 
1 and is shown schematically in Fig. 1. The gas state proper- 
ties and the component flowfield models are discussed in detail 
in the following subsections. The results of applying the models 
to three Jovian entry conditions are presented in Section 3. 

2.1 Thermodynamics 

Equilibrium thermodynamic properties were used throughout 
this study. The curve fit coefficients for heat capacity, en- 
thalpy and entropy were obtained from Moss (Private communica- 
tions from James Moss, NASA Langley Research Center, Hampton, 
Virginia) for use in the chemical equilibrium composition (CEC) 
program (Ref. 2). Two ranges of thermodynamic property curve 
fits were used (1000 K to 6000 K and 6000 K to 16000 K). The CEC 
program determines the equilibrium composition of a gas mixture 
by minimization of Gibbs free energy. 

The three flight cases analyzed during this study are speci- 
fied in Table 1 . Case 4 and 5 are documented in Ref. 3 and Case 
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3b was obtained from Moss (Private Communications from James 
Moss, NASA Langley, Hampton, Virginia). Computer printout from 
the HYVIS program (Ref. 3) were obtained for all three cases. 
Thermodynamic property output from the Langley HYVIS program were 
compared with the use of the CEC program in this study. These 
results are shown in Tables 2a, b and c for Cases 3b, 4 and 5 
respectively. The temperature, pressure and elemental composi- 
tion were specified and all other properties were computed by the 
CEC program. The results are for three shock layer conditions 
and are in excellent agreement for all variables. 

The CEC program was used to calculate thermodynamic property 
information for the flowfield math models used in the wake. To 
accomplish this several options were used: 

RKT - Constant entropy expansion from specified total 
conditions 

T,P - Specified temperature and pressure 

H,P - Specified enthalpy and pressure 

Throughout the analysis six species were considered in the invis- 
cid flowfield and nineteen species were considered in the viscous 
f lowf ields . 

2.2 Shoulder Expansion 

In order to determine the initial conditions for the near 
wake viscous shear layer, an isentropic stream tube approximation 
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j 

I 

I 

was used to calculate the boundary layer flow over the probe 
shoulder. The computational area of interest is shown in Fig. 2 

i 

, and is confined by the body and outer edge of the viscous region 

over the probe shoulder. 


The geometric relations used to determine the areas before 

and after expansion are shown in Fig. 3. HYVIS output for 

all points in the boundary layer are used along the 

Y-coor di na te . These points are considered as stream tubes and 

expanded to locations along the Z-coordinate . The angle, 0 , is 

w 

the wake angle determined by the method of characteristics (MOC) 
solution just downstream of the shoulder. To expand the flow ar- 
ound the shoulder the following relations are used: 

Known Upstream Known Downstream 

Conditions Conditions 


Ui 

P. 


from HYVIS 


Z = 0 

- from MOC 


Using continuity we may write: 

^ 1+1 - * 1+1 

To determine velocity the adiabatic perfect 


_L in. n -x(i) 

2 2 Y-l\p/ 


( 1 ) 

gas equation is used. 

( 2 ) 


Where the isentropic perfect gas relation for pressure and densi- 
ty is required. 

P/P » (P/P)Y 
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( 3 ) 



The area relations from Fig. 3 are used to locate the stream- 


lines 

-= TT(2RB+(Y^^.l+Yi)cosec)[Yi+i-Y^] 
^i+1 “ TT(2RB+ri+i+ri)[?i+i-ri]/cos0„ 

Let R = Rg + f and solve for rj[+i 
Xi+i = TT(R|4.i-Ri)/cos6„ 

^i+1 “ (cos6„ ~ % 


(4) 

(5) 

( 6 ) 
(7) 


The key point in the expansion method is that the boundary 
layer must be expanded to the pressure defined by the MOC calcu- 
lation. The pressure gradient predicted by the MOC near the wall 
is quite large as shown in Fig. 4 for the cases analyzed. These 
profiles were imposed on the viscous layer and yield realistic 
viscous edge conditions after expansion as shown in Fig. 5. The 
theoretical edge locations are in reasonable agreement with the 
experimental data from Park, (unpublished data from NASA AMES) . 
However, the ablation rate for the ballistic range test of Park 
are unknown and initial flow conditions are different. 


2.3 Inviscid Near and Far Wakes 

The axisymmetric methods of characteristic (MOC) code of 
Ref. 4 as modified in Ref. 1 was used to calculate the flow 
properties of all the inviscid regions shown in Fig. 1. The 
code operates with entropy and Mach number as independent thermo- 
dynamic variables for a constant total enthalpy. The code for 
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the Jovian probe problem uses thirty velocities at each of the 
thirty entropy cuts. The entropy and velocity values were se- 
lected for each case considered. The MOC code was modified to 
provide a summary print of flow properties along the shock, the 
body surface and along any free boundary. 

The probe’s forebody was modeled as an equivalent cone. The 

forebody angle , 9 , was altered to obtain an approximate match of 

c 

post shock value with those of the HYVIS output. The forebody 

was followed by a short cylindrical section (probe shoulder). 

The cylindrical section was followed by a constant pressure free 

boundary. The base pressure correlation shown in Fig. 6 from 

Ref. 1 was used to establish the pressure level. The pressure 

values for the cases previously studied (Ref. 1) along with the 

values for the current cases are shown in Fig. 6 . The lower 

boundary correlation was used to get the neck location within the 

neck location bounds developed in Ref. 1. These boundaries are: 
Neck Location 

X/R - 2.8 ± 0.4 ( 8 ) 

r/R » 0.4 ± 0.2 

Downstream of Neck 

r/R « b (X/R)® ( 9 ) 

m «• 0.25 ± 0.05 

The preceeding relations were based on thin boundary layer wind 
tunnel data. In the current problem a blown boundary layer ex- 
ists which significantly alters the effective body geometry. 
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Consequently, an effective body radius which extends to the outer 
streamline of the viscous layer, see Fig. 5, was used to replace 
the body radius, R, in Eq . 8 and 9* This permitted obtaining 
consistent neck locations from Eq . 8 and 9 and from the MOC 
using the pressure correlation. 

The far field inviscid wake started with initial conditions 
from the near wake calculation along the vertical dashed line 
shown in Fig. 1. The MOC treated the far viscous wake as a 
solid boundary defined by Eq. 9* The wake shock was calculated 
as part of the overall MOC solution. Calculations were made from 
the wake neck to X/R = ?• 

The results from the MOC calculation are put on tape and are 
arranged along characteristic lines in cylindrical coordinates. 
A program was written to arrange the data in uniform radial in- 
crements at specified axial values. This program called JUGGLE 
also permits translation of coordinates from those used by the 
MOC to those desired for final data usage. The JUGGLE program 
was used on both the near and far inviscid wakes. 

2.4 Viscous Near and Far Wakes 

The near and far viscous wakes were calculated using a modi- 
fied form of the computer program documented in Ref. 6. The 
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basic features of the program and its current application are 
described below. 


The conservation equations for two dimensional or axisymme- 
tric viscous mixing flow were written in terras of stream func- 
tions. An explicit finite diffrence equation was constructed for 
the axial momentum equation to evaluate velocity. The axial mom- 
entum equation is: 


8u _L ^ 1 ^ / 9u \ 

9x “ ” pu dx ^ 3«(;! j (10) 

where a =• ypuy^*^/^;® (11) 

and ■ puy" (12) 

The term n is 0 for two-dimensional flow and 1 for axisym- 
metric flow. The resultant velocity values are applied in a sim- 
ilarity procedure to determine the degree of mixing. The degree 
of mixing is calculated using 



(13) 


In the current application, F represents the fraction of abla- 
tion products mixed with the external H 2 and He species and the 
nondimensional static enthalpy level. The similarity approxima- 
tion can be assessed in terms of the HYVIS output at aft end of 
the forebody conical section (station 17). Results for the three 
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cases analyzed are shown in Fig. 7 and indicate that the simi- 
larity approximation provides a reasonable set of initial condi- 
tions for the viscous near wake. 


The temperature of the mixed composition is then determined 
under the restraint of chemical equilibrium conditions. The vel- 
ocity, temperature and other properties are located in the physi- 
cal plane by an inverse transformation from the stream function 
plane . 


The eddy viscosity model used in this study is a 
Prandtl-like model for compressible flow: 


initial 

900 


(pU) - (pU) . 1+10 

max mxn ' 


-4 


(14) 


Other models were examined but were found to yield more rapid 
spreading than the observed viscous wake boundaries would indi- 
cate is realistic or no significant difference from the preceed- 
ing model. 


The chemical equilibrium shear layer (CESL) program was used 
to calculate the near wake and far wake viscous areas. Both near 
and far wake calculations required special handling as described 
below . 
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The initial conditions for the near viscous wake was provid- 
ed by the shoulder expansion code described in Subsection 2.2. 
The key concept used in obtaining a realistic expansion of the 
viscous layer was to superimpose the pressures determined from 
the MOC calculation. This same concept was extended to the outer 
portion of the viscous near wake as schematically shown in Fig. 
2. An example of the pressure profiles derived from MOC calcula- 
tions for Case 3b is shown in Fig. 8. The outer edge of the 
shear layer is also shown in Fig. 8 and illustrates the extent 

to which the viscous layer is modified by the imposed MOC pres- 
sure field. 

The second important concept used in developing the near 
viscous wake was specifying the outer wake streamline location. 
This method eliminates the need to transform two dimensional so- 
lutions as done by Balakrishnan and Chu (Ref. 7). it properly 
accounts for the axisymmetric convergence of the flow. Further, 
the recompression pressure in the near wake is calculated as part 
of the solution. This type of approach is necessary since the 
shear layer is not thin compared to any dimension of the body or 
wake. In order to specify its physical location in the inner re- 
pressure must rise to accommodate the area being restrict- 
ed. Thus, an axially varying constant-radial-pressure distribu- 
tion is calculated. A forward integration from shoulder to wake 
neck region is made and a neck pressure is calculated. This 
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pressure is compared with the MOC predicted pressure. If they 
are not equal, the neck radius is changed and a new integration 
is begun. 

Experimental and theoretical pressure distributions for tur- 
bulent wakes were examined to verify the expected trends. Figure 
9 from Ref. 8 shows the pressure distribution for a 
two-dimensional case. Note that there is no constant pressure 
region in the data. Figures 10 and 11 from Ref. 9 show pressure 
distributions for axisymmetric turbulent wakes. Figure 10 is a 
comparison of experimental and theoretical pressures along the 
axis of symmetry. Both theory and data are for a cylindrical 
forebody. Peak pressure on these graphs and the wake neck occur 
at the same axial location. The theoretical results given in 
Fig. 11 indicate that a thick initial boundary layer 
thickness, 5j, yields a very short constant pressure region. 

In the present problem we have both a large 6i and a high 
expansion angle. The expansion angle is twice that of Fig. 9. 
Thus, no significant constant pressure region is expected. The 
results shown in Section 3 from the current model appear consis- 
tent with the available data and expected trends. 

The far wake calculations were made with an axial pressure 
and corresponding edge condition variation predicted by the MOC 
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results. 


These calculations were made with axisymmetric equa- 
tions without boundary location constraints. Initial conditions 
for the far wake were taken directly from the near wake neck lo- 
cation results. 

2.5 Recirculation Zone 

Defining the recirculation zone properties for the Jovian 
probe presents a rather unique base flow problem. Several as- 
pects of this problem make it unique. First, very little experi- 
mental work has been done on turbulent wakes of large angle 
cones. Secondly, most experimental and theoretical wake studies 
have been done on thin boundary layers and not blown boundary la- 
yers as exhibited by the Jovian probe. Thirdly, few studies have 
addressed determining the recirculation zone enthalpy level. 
Fourth, high Mach number data for turbulent wakes are scarce. 
Finally, finding information on axisymmetric turbulent wakes 
which contain more than one of these aspects is rare. 

In order to determine the enthalpy and concentration levels 
in the recirculation zone, a physical and resulting math model 
must be adopted. Figure 12 provides schematics for five physical 
models examined during the course of this study. Model 1 has a 
closed dividing streamline and is the model usually found in the 
literature. The model assumes steady state flow with the mass in 
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the recirculation zone constant. Problems were found with using 
Model 1 for the current work. First, since the mass in the re- 
circulation zone cannot be exchanged with the external shear 
layer, cannot be determined. Secondly, if a steady state energy 
balance is made at the boundaries of the zone in Model 1 , unreal- 
istic results are obtained. If the body is assumed adiabatic, 
the recirculation zone total enthalpy is the freestream total en- 
thalpy. This implies that the forebody massive ablation provides 
no protection to the base. If the body is nonadiabatic then the 
energy transferred across the dividing streamline is transferred 
to the body. Using this assumption the enthalpy level cannot be 
determined. Based on the foregoing reasons Model 1 was rejected. 

Models 2 and 3 are similar in that one end of the dividing 
streamline is open to mass flow and steady flow is assumed. 
Since Model 2 would lead to an accumulation of mass and Model 3 a 
depletion of mass within the recirculation zone, both models were 
rejected. 

Model 4 assumes steady flow and an open recirculation zone 
in which some mass from the boundary layer is captured and the 
same amount is allowed to exit at the wake neck. The recircula- 
tion zone is like a stirred reactor with a constant mass input 
and output. The mechanisms which would allow both a capturing 
and escape of mass could not be identified from the literature. 
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If this process is real, a capturing model for flow over the 
shoulder would require development. 

The literature on separated and wake flows by Chilcott (Ref. 
10) and Charwat et. al. (Ref. 11) suggest that a pulsating 
mass exchange process exists and the degree to which this behavi- 
or dominates is geometry dependent. Using this concept, Model 5 
was developed. The dividing streamline is closed at the shoulder 
and periodically open at the wake neck. The opening process 
loses mass and the closing process pumps mass into the recircula- 
tion zone. Charwat et. al. (Ref. 11) points out that even for 
cavity type flows neither the "conduction” model, corresponding 
to Model 1 , nor the mass exchange model are sufficiently com- 
plete. In the current problem, effects of Model 1, 4 and 5 are 
thought to contribute to the recirculation zone properties. The 
conduction across the dividing streamline, and the mass exchange 
mechanisms of Models 4 and 5 determine the recirculation zone en- 
thalpy. The mass exchange mechanisms of Models 4 and/or 5 deter- 
mine the elemental composition within the recirculation zone. 

Since development of a steady state model based on global 
energy and mass balances do not appear feasible and correct, an 
empirical approach was taken. Experimental data exist for the 
recirculation zone temperatures of nonablating bodies. Huber and 
Hunt (Ref. 12) compiled data to compare with the Apollo shape. 
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These data are shown in Fig. 13 where temperature was replaced 

by enthalpy. Turbulent flow data corresponding to Re >10® for 

the Apollo shape was selected to best represent the Jovian probe. 

Wall enthalpy does not Bsignif icantly influence the recirculation 

zone enthalpy since H /H„=0.030 to 0.039 for the cases analyzed. 

w T 

A value of (H, -H ) = 0.5 was selected as a conservative 

b w T w 

estimate for the design application herein. An improved estimate 
and model awaits test data for hypersonic blunt cones with mas- 
sive forebody blowing. 
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Section 3 
RESULTS 


This section presents observations about the three wake 
flowfields calculated and presents the format in which the flow- 
fields were recorded for radiation analysis work. 

3.1 Flowfield Distributions 

The analytical methods described in Section 2.0 were used to 
calculate three flowfields corresponding to three possible Jovian 
probe entry conditions. The specifying conditions for the three 
cases are given in Table 1. Case 3b corresponds to peak forebody 
heating for a cool atmosphere model with a steep entry angle. 
Cases 4 and 5 are for a nominal atmosphere model and nominal 
entry angle (Ref. 3). Case 4 corresponds to peak forebody heat- 
ing in the nominal entry trajectory. 

The starting location for the current flowfields is where 
the HYVIS program terminates at the end of the Jovian probe coni- 
cal section (Station 17). It is therefore appropriate to compare 
the MOC solution on the conical section with the HYVIS results. 
MOC results are for a sharp cone whereas the HYVIS results are 
for a sphere-cone body. Post shock results are given in Table 3 
and are found to be in good agreement. The body angle used in 


17 


the MOC calculation was larger than the actual body angle to ac- 
count for effective body displacement due to massive blowing. 

The post shock conditions shown in Table 3 are indicative of 
the relative temperature and pressure levels in the respective 
wakes. This point will become apparent as the results given in 
this section are reviewed. The post shock conditions for the 
forebody shock past the probe’s shoulder are shown in Fig. 14. 
Case 3b exhibits the highest temperatures whereas Case 5 exhibits 
the lowest temperatures. This same trend persists throughout the 
near and far wake inviscid flowfields as summzrized in Fig. 15. 
Post wake shock conditions may be significant for Case 3b but ap- 
pear somewhat benign for Case 5. 


Streamline pressures and temperatures for the viscous near 
and far wakes are shown in Figs. 16 to 21. The external and 
internal viscous streamline pressures, shown in Figs. 16, 18 and 
20, converge at X/R ~ 0.5. Between X/R = 0 and 0.5 large radial 
pressure gradients exist in the shear layer as shown in Fig. 8. 
This high pressure and pressure gradient region produces the high 
temperature levels shown in Figs. 17, 19, and 21. The models 
used in the previous work (Ref. 1) did not account for this 
pressure effect and thus did not identify this high temperature 
region. The rise in temperature at the neck location is a result 
of the wake shock. Figures 17, 19, and 21 also show the recircu- 
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lation region temperature levels. The slight rise with increas- 
ing X/R is due to the axial pressure rise in the recirculation 
zone which matches the inner streamline pressure. 

Near wake radial profiles for the peak forebody heating con- 
ditions of two trajectories are shown in Fig. 22. Pressure lev- 
els and temperature levels in the recirculation region, viscous 
layer and inviscid region are higher for Case 3b than for Case 4. 
Mass fraction levels of carbon in the recirculation and viscous 
layers are comparable for the two cases. 

The nondimensional pressure distributions in the near w,ake 
calculated using the model described in Section 2.3 are shown in 
Fig. 23- The effect of massive blowing increases the boundary 
layer thickness and nearly eliminates the plateau pressure region 
in the near wake. This trend has been illustrated for less sev- 
ere cases by Ref. 9 as shown in Fig. 11. The base pressure is 
increased but only slightly by increased ablation. The calculat- 
ed neck pressure for the three cases is between 2.5 and 3.3 times 
the base pressure. This agrees with the data and theory of Fig. 
10 and data and theoretical results presented by Wong and Chow in 
Ref. 21. However, the current results are significantly differ- 
ent from those presented by Park (Ref. 22). 

Park in Ref. 22 assumed a base flow model like Model 1 of 
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Fig. 12 to establish the neck pressure. A dividing streamline 
velocity of 50 percent of the external velocity was assumed to 
stagnate and close the recirculation zone. This leads to the 
high neck pressures presented in Ref. 22. No explanation of how 
the recirculation mass fraction or enthalpy level can be estimat- 
ed using the closed recirculation zone concept was given. Park 
estimated the neck to stagnation pressure ratio = 0.15 for 
Case 4 of this report. This may be compared to the current re- 
sults in Table 4. By using base heating rate correlation of Ref. 
22 with the current neck pressures, the base radiative heating 
was estimated to be 0.925 MW/m^ compared to the value of 19-4 
MW/m^from Ref. 22. The actual value for Case 4 from the current 
analysis must await the detailed radiative calcualtions for the 
flowfield provided by this work. However, by using the correla- 
tion of Park the current results would appear to yield heating 
rates lower by a factor of 20 than results of Park's and higher 
by a factor of about 3 than the results from Brant and Nestler 
(Ref. 5). 

The main uncertainty and differences between different ana- 
lysis is the manner in which the recirculation region is closed. 
This produces different models for the recirculation zone species 
and enthalpy level as well as different wake neck pressure lev- 
els. All of these factors significantly affect the base heating 
environment. This uncertainty will most likely persist until an 
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experimental data base is available for blunt planetary probe 
class bodies and entry conditions. 

3.2 Flowfield Format 

The flowfield results for all three cases were delivered to 
NASA-Langley on cards. The flowfield properties were provided in 
nondimensional cylindrical coordinates starting at the aft end of 
the probes shoulder and going to X/R = 7 « 0 . The card format is 
defined in Table 5 . The flowfield for each case consists of ap- 
proximately 1300 points where temperature, pressure and elemental 
composition is defined. 
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Section 4 

CONCLUSIONS AND RECOMMENDATIONS 


Based on the analysis and results of this study, the follow- 
ing conclusions are drawn. 

(1) Case 3b exhibited the highest wake temperatures 
and Case 5 the lowest wake temperatures. 

(2) The neck pressure was calculated to be 
approximately 3 times the base pressure. 

(3) The current pressure model for the outer 
portion of the shear layer produces a high 
temperature region in the shear layer for 
X/R 0.5. 

(4) The recirculation zone enthalpy can be only 
approximated using current models and 
experimental data. 

(5) Large differences exist in the neck pressures 
used by different options to model the Jovian 
probe wake. 


As a consequence of the information obtained in this study, 
the following recommendation is made. 

An experimental program, for general planetary probe 

application, should be undertaken to establish: 

(a) The pressure distribution in the near wake both 
with and without massive blowing. 

(b) The elemental composition in the recirculation 
zone as a function of blowing rate. 

(c) The recirculation zone enthalpy as a function of 
blowing rate and mass injection enthalpy ratio to 
total enthalpy. 
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Table 1 


SPECIFYING CONDITIONS 


Variable 

Case 3b 

Case 4 

Case 5 

Time (sec) 

85.25 

111.3 

113.5 

Altitude (km) 

93.346 

126.05 

115.31 

Atmosphere Mass Fractions 




H2 

.641 

.803 

.803 

He 

.359 

.197 

.197 

Poo (atm) 

.004429 

.002422 

.003795 

Too (°K) 

132.07 

152.00 

148.0 

Moo 

46.975 

43.009 

38.45 

Voo (m/sec) 

37547 

39284 

34666 
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Variable 


P atm 
T K 

H cal/gm 
M Mol wt 


Table 2a 
Case 3b 

COMPARISON OF EQUILIBRIUM THERMODYNAMIC CALCULATIONS 


Condition 3 


12.651 

17125 


1.118 

26.2726 

1.2215 

12474.3 

.52184 

.00003 

.18890 

.10016 

.00009 

.18899 


12.651 

17125 


166386. 166322. 


1.118 


.52199 

.00003 

.18890 

.09999 

.00009 

.18899 


Condition 2 

REMTECH 

Langley 

7.1249 

7.1249 

12999 

12999 

94742 

94716 

1.307 

1.307 

12.7979 


1.2480 


10158.7 


.77931 

.77949 

.00010 

.00009 

.05168 

.05169 

.11723 

.11704 

.05168 

.05169 



7.436 

3978 

5417.9 

30.824 

2.1459 

1.1076 

1090.1 

.10850 

.03768 


7.436 

3978 

5423.9 

30.787 



.10904 

.03768 


.01757 .01752 



Condition 

1 - Post Shock (Stagnation Point) 

2 - Post Shock (Station 17) I .. .. j r r u j i 

3 - Wall (Station 17) 1 At aft end of forebody comcal section 
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Table 2b 
Case 4 

COMPARISON OF EQUILIBRIUM THERMODYNAMIC CALCULATIONS 



Condition 1 

Condition 2 

Condition 3 

Variable 

REMTECH 

Langley 

REMTECH 

Langley 

REMTECH 

Langley 

P atm 

6.008 


3.356 

3.356 

3.469 

3.469 

T “K 

15742 

15742 

11822 

11822 

3878 

3878 

H cal/gm 

182219.8 


102858.3 

102838.9 

5510.9 

5506.5 

M Mol wt 

0.972 

0.972 

1.137 

1.137 

29.986 

30.009 

Cp cal/gmK 

32.7970 

— 

13.3459 

— 

2.4762 

— 

Y 

1.2105 

— 

1.2560 

— 

1.1071 

— 

a m/sec 
Mole Fractions 

12767.6 

— 

10421.4 

— 

1091.0 

— 

H 

0.59623 

.59845 

0.86674 

.86683 

0.13143 


Hz 

0.00002 


0.00009 

.00009 

0.03700 


H+ 

0.17794 

.17792 . 

0.03860 

.03861 

— 


He 

0.04783 

.04776 ■ 

0.05596 

.05587 

0.00742 

.00743 

He+ 

0.00002 

.00002 

— 

— 

— 

— 

e- 

0.17796 

.17794 

0.03860 

.03861 

— 

— 

C 

— 

— 

— 

— - . 

0.00968 

.00960 

Cz 

— 

— 

— 

— 

0.01156 

.01150 ^ 

C 3 

— 

— 


— 

0.20350 

.20291 

CjH 


— 

— 

— 

0.12161 

.12157 


— 

— 

— 

— 

0.01922 

.01930 

C 3 H 

— 

— 

— 

— 

0.13837 

.13863 

C,H 

— 

— 

— 

— 

0.20833 

.20905 

CO 





0.11187 

.11199 


Condition 


1 - Post Shock (Stagnation Point) 

2 - Post Shock (station 17) 

3 - Wall (Station 17) 
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Table 2c 


Case 5 

COMPARISON OF EQUILIBRIUM THERMODYNAMIC CALCULATIONS 



Condition 1 

Condition 2 

Condition 3 

Variable 

REMTECH 

Langley 

REMTECH 

Langley 

REMTECH 

Langley 

P atm 

7.447 

7.447 

4.105 

4.105 

4.228 

4.228 

T 

14475 

14475 

8839 

8839 

3765 

3765 

H cal/gm 

141561.3 

141531.1 

77711.9 

77700.0 

5498.5 

5503.9 

M Mol wt 
Cp cal/gmK 

Y 

1.059 

22.929 

1.2187 

1.059 

1.180 

5.2033 

1.5199 

1.180 

23.600 

2.4672 

1.1163 

23.579 







a m/sec 
Mole Fractions 

11766.9 

9729.8 

1216.9 




H 

.73960 

.73967 

.93601 

.93611 

.18329 

.18400 

Ha 

.00006 

.00006 

.00061 

.00061 

.13453 

.13427 

H+ 

.10410 

.05214 

.10411 

.05205 

.00265 

.05809 

.00265 

.05799 



He 

He+ 

e- 

.03299 

.03291 





.10410 

.10411 

.00265 

.00265 



C 

.00396 

.00399 

Ca 





.00412 

.00416 

C3 





.07108 

.07146 

CjH 





.12677 

.12690 

^2^a 





.05255 

.05229 

C 3 H 




1 

.13260 

.13242 

C.,H 





.17211 

.17165 

CO 





.08599 

.08593 


Condition 1 - Post Shock (Stagnation Point) 

2 - Post Shock (Station 17) 

3 - Wall (Station 17) 
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Table 3 

COMPARISON OF CONE SECTION POST-SHOCK CONDITIONS 




Case 3b 

Case 4 

Case 5 

Variable 

Units 

REMTECH 

LANGLEY 
Sta. 17 

REMTECH 

LANGLEY 
Sta. 17 

REMTECH 

LANGLEY 
Sta. 17 

T 

K 

13,044 

12,999 

11,908 

11,822 

8994 

8839 

P 

atm 

7.241 

7.125 

3.432 

3.356 

4.213 

4.105 

Avg. Mol . 
Weight 

— 

1.306 

1.3074 

1.1345 

1.137 

1.179 

1.180 

es-6c 

(Blowing) 

deg. 

3.91 

4.18 

3.733 

3.88 

3.47 

2.902 

6c 

deg. 

46 

45 

46 

45 

45.4 

45 
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Table 4 


BASE PRESSURES SUMMARY AND BASE RADIATIVE HEATING ESTIMATES 


Case 

Ps 

(atm) 

— 

Pn/Ps 

PN/Pb 

V'ls 

(1) 

ti) 

(MW/m^) 

IHi 

Ha 

■oPI 

.2721 

.0215 

3.32 


451.9 


HI 

BEH 

.0964 

.0160 

2.52 

.00488 

189.5 

0.925 


7.447 

.1278 

.0168 

2.59 


166.3 

0.605 


Notes (1) = 2.4 (Pj^|/Pg)^® from Park (Ref. 22) 

(2) q = stagnation coupled ablation radiative heating 

from HYVIS 

(3) qjj = base radiative estimate using (1) and (2) 


Table 5 


CARD FORMAT 


Col umn 

Format 

Variable 

1 - 10 

Flo. 6 

X/R, nondimensional axial distance 

11 - 20 

FI 0.6 

r/R, nondimensional radial distance 

21 - 30 

Flo. 4 

Static temperature (®K) 

31 - 40 


Static pressure (atm) 

41 - 50 

F10.7 

Mass fraction of hydrogen 

51 - 60 

F10.7 

Mass fraction of helium 

61 - 70 

F10.7 

Mass fraction of carbon 

71 - 80 

F10.7 

Mass fraction of oxygen 
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^ 

Cone Frustum Areas 

Ai+i = IT ((r^+i+Re) + (ri+RB))[Y^+i- Y^] 

Ai+i = TT ((r-j+l+Re) * (?‘i'‘’RB))CZi+l“Zi] 

r = Y COS0J, 

r = Z cos6y, 

Therefore 

Ai*+i = TT (2RB+(Yi+i*Yi)cos6c)[Yi+i-Y^] 

A^+i = TT (2RB+ri+i+r^)[r^+l-r^]/cos0w 

Fig. 3 Shoulder Expansion Geometry 
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Pressure (Ib/ft^) 



r/R r/R r/R 

Case 3b Case 4 Case 5 

Fig. 4 Shoulder Aft Position Pressure Profiles 



Fig. 5 Boundary Edge Geometries Over the Probes Shoulder 
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Fig. 6 Turbulent Base Pressure Correlation for Planetary 
Entry Systems (From Brant and Nestler, Ref. 5) 
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Pressure (Ib^/ft^) 
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O static pressure probe 



Fig. 9 Two-Dimensional Turbulent Wake Spliter Plate Pressure Distribution 
(Moo ■ 20.4, Rboo = 3.4x10®, 0V7=22.5“) [From Wagner (Ref. 8 )] 
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omo] 


imi 


Mp2.0 
ei/R = 0.0061 
EXPERIMENT 

(Reid & Hastings) 


M]=1.92 

Gi/R=0.02405 

EXPERIMENT 

(Sirieix, Delery & Monmerie) 
o o 


Mpl.85 

e^/R=0.05575 

EXPERIMENT 

(Baebinarayanan) 


10 Comparison of Experimental and Theoretical Turbulent 
Axi symmetric Wake Centerline Pressure Distributions 
(From Mehta and Strahle, Ref. 9) 




Fig. n Comparison of Theoretical Turbulent Axisymmetric Wake 
Centerline Distributions for Different Initial Boundary 
Layer Thicknesses (From Mehta and Strahle, Ref. 9) 
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Rear Stag. Point 

imrir cones, wedges 0.2 13,14,15,16 

circular cylinder 0.2 17 

Recirculation Region 
- wedges 0.125 18 

Gemini , a=-15* 0.40 19 

□ cone 0.58 20 

Apollo, a=20« 0.40 12 



I i— _1_J i 1 I I I— I 

0.0 .1 .2 .3 .4 .5 .6 .7 .8 


(Hb-Hw)/(Hx-Hw) 


Fig. 13 Comparison of Reynolds Number Dependence on Near-Wake 
Enthalpy (From Ref. 12) 
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r/R 


Case 

M 




3b 

cx>= 43.95 



IN ITALICS ARE P/Pco) 


3280 K 
* 95. S 


3680 K 
’ 30. S 


3710 K 



^ 3130K 
^«2840K’ 83. S 
188 


3400K 

29.2 


3550K 


3465K 

*16.3 



2900K 
• 77.1 
2615K 
149 

. 3560K 

22.0 


3380K 3640K 

12.7 86.4 




3675K 

29.0 


3260K 


3590K 


3590K 

23.4 


13.7 

25.0 


t - - > - 1 » * ^ > ■ 

3 4 5 6 7 

X/R 


15 Comparison of Inviscid Wakes 





Pressure (atm.) 



X/R 

Fig. 16 Viscous Wake Streamline Pressures for Case 3b 




X/R 

Fig. 17 Viscous Wake Streamlines and Recirculation Zone Temperatures for Case 3b 





X/R 

Fig. 18 Viscous Wake Streamline Pressures for Case 4 
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Fig. 19 Viscous Wake Streamlines and Recirculation Zone Temperatures for Case 4 
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:;j Inner Viscous Streamline : i 


3 4 5 6 7 

X/R 


Fig. 21 Viscous Wake Streamlines and Recirculation Zone Temperatures for Case 5 
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Fig. 23 Near Wake No 
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